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Abstract 

Charmed meson charmonium spectra are studied with improved quark actions on 
anisotropic lattices. We measured the pseudo-scalar and vector meson dispersion 
relations for 4 lowest lattice momentum modes with quark mass values ranging 
from the strange quark to charm quark with 3 different values of gauge coupling (3 
and 4 different values of bare speed of light v. With the bare speed of light parameter 
v tuned in a mass-dependent way, we study the mass spectra of D, D s , r/ c , D*, D* 
and J/ijj mesons. The results extrapolated to the continuum limit are compared 
with the experiment and qualitative agreement is found. 
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1 Introduction 

In recent years, anisotropic lattices and improved lattice actions have been 
used extensively in lattice QCD calculations. In this letter, we present our 
preliminary numerical study on the spectrum of charmed mesons and char- 
moniums. The gauge action employed in this work is the tadpole improved 
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gluonic action on asymmetric lattices as in Ref. [2]: Using such a gauge ac- 
tion, glueball and light hadron spectrum has been studied within the quenched 
approximation [1,2,3,4,5,6,7]. 

It has been suggested that relativistic heavy quarks can also be treated with 
the help of anisotropic lattices (the Fermi lab approach), with improvements 
[8,9,10,11,12]. Using various versions of the quark actions, charmed meson 
spectrum, charmonium spectrum, charmed baryons have been studied on the 
lattice [13,14,15,16,17,18]. However, in order to take full advantage of the 
improved quark action on anisotropic lattices, some parameters in the action 
have to be tuned, either perturbatively or non-perturbatively, in order to gain 
as much improvement as possible. Some numerical studies of these parameters 
have already appeared in the literature [19,12,20,21,22]. In particular, it turns 
out that the bare speed of light parameter v has to be tuned non-perturbatively 
in a quark mass- dependent way as is shown in Ref. [23]. With the parameter 
v tuned as such, we can study quenched meson spectrum with quark mass 
values ranging from the strange to the charm, using the same quark action. 

The anisotropic quark actions used in our study is: Sf = J2 xy i^x-M-xy^y with 
the fermion matrix given by [23]: 



1/ (2K max ) + Pt^CToiFoi + P s [(JyiTy1 + 0-23^23 + 0"31-^3l) 
i=l 

\(1 - 7»)Uv(x)6 x+ w + (1 + i^Ufa - riSz-w] . (I) 



The coefficients in the matrix are given explicitly by: 
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In this letter, we present our preliminary quenched results for the meson spec- 
tra using the above improved quark action on anisotropic lattices. 



2 Simulation Results 



The whole procedure of the calculation involves three major steps: In the first 
step, we determine the optimal value of the bare speed of light parameter v, 
using energy-momentum dispersion relations of pseudo-scalar mesons; In the 
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second step, the masses of the mesons that we are interested in are extrap- 
olated/interpolated to the optimal value of is; In the third step, we perform 
chiral and continuum extrapolations as usual. 

The lattices in this study are of size 6 x 9 x 12 x 50 for (3 = 2.4,2.6 and 
8 x 12 x 16 x 50 for (3 = 2.8. The spatial lattice spacing a s roughly range from 
0.13/m and 0.23 fm in physical units. The aspect ratio is £ = a s /a t ~ £o — 5 
for all lattices. For gauge field configurations at a given value of (3, 4 different 
values of the bare speed of light parameter is and 12 different values of the 
hopping parameter k are studied. The range of k roughly covers the quark 
mass from the strange to the charm. About 300 configurations are used for 
each parameter set in this study. 

The energy values of meson with definite three-momentum k (including zero- 
momentum) is obtained from their respective correlation functions by finding 
the plateaus in their effective mass plots. The errors for the data points in the 
effective mass plots are analyzed using the standard jack-knife method. We 
obtain from the pseudo-scalar and vector channel the corresponding energy 
Eps/v( m i, m 2, v\, is 2 , k), where rrt\ and m<i are the bare quark mass parameters 
defined by: m = l/(2/c) — l/[2/t cr (z/)]. The value of l/[2/c cr (i/)] is obtained by 
the condition of vanishing pseudo-scalar meson mass for each is. Although for 
different values of is, the set of values for m turns out to be slightly different, 
we always interpolate/extrapolate them to the same set of values [23]. 

To obtain the optimal value of is, the values of Ep s (m, m, is, is, k) are fitted 
linearly against the three-momentum squared k 2 in the low-momentum region: 

Eps( m ' m i v i v i k) = Mp S (m, m, is, is) + Z PS {rn, m, is, is)k 2 . (3) 

The condition: Z PS [m, m, is opt {m) , v pt( m )] — 1 then yields the optimal value 
of v as a function of the bare quark mass parameter: v opt {rn). In practice, we 
use the values of Zps(m, m, is, is) at different is to perform a linear extrapola- 
tion/interpolation in is for each value of m. 

We now proceed to the second step of our study, namely to extrapolate/interpolate 
the meson mass values to the optimal values of is obtained from the pseudo- 
scalar dispersion relations. This is achieved by a linear extrapolation/interpolation 
in is, or equivalently in the quantity Z PS {m, m, is, is). 

Similar procedure discussed above can be repeated for mesons consisting of a 
quark and an anti-quark with different mass. We interpolate/extrapolate the 
meson masses to the optimal bare speed of light parameter. The outcome of 
this procedure are the quantities of Mp s , v (mi,m2,is pt(mi),is opt ('m2))- Since 
these quantities only depend on the bare quark mass parameter pairs (mi, 1712), 
we denote them by Mp S ^(mi,m 2 ). 
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In the third step of the analysis, we perform the chiral extrapolations. First, 
we determine the physical strange and charm bare quark mass parameters: 
m {phy) an( j m (phy) _ j n fais work, we use the mass of the vector meson with 
the same flavor of quarks to determine m^ hy ^ and m < f hy \ The results of 
My(m, m, f op t(m), f op t(m)) are fitted quadratically versus the quark mass pa- 
rameter: 

afMy[m, m, ^ op t(m), is opt (m)] = A + Bin + Cm 2 . (4) 

We perform two such quadratic fits, one in the low quark mass region, the 
other in heavy quark mass region. The scale of the lattice is set by r = 
0.5fm (the so-called Sommer scale). Then the value of m!f hyS) and m t f hy " > are 
obtained by requiring the vector meson mass obtained from the simulation 
to exactly reproduce the mass of the physical $ and the physical J/ip meson 
mass. Fig. 1 shows the quadratic extrapolations of M v (m, m, u opt (m), ^ op t(m)) 
versus the bare quark mass m at (5 = 2.6. The physical $ and J/psi mass 
are also lined out by horizontal lines in the figure. The intersection point of 
the lower/higher horizontal line representing a 2 M| / af M 2 ^ and the quadratic 
fitting curve yields the estimate for the physical strange quark mass parameter 
m (phy) anc [ m (p h y)^ respectively. 

Having determined the values of m < f hy "> and then make quadratic ex- 

trapolations/interpolations for all the meson mass squared (Mp S ^ v (m 1 ,m2)) 
versus the quark mass parameters and extrapolate/interpolate them to their 
physical values rnfj 1 ^ , depending on which type of meson is being considered. 
The final results are the mass values of various charmed mesons and char- 
monium states at different (3. These results are listed in Table tab:spectra. 2 

Table 1 



Extracted physical mass (in unit of MeV) of vector and pseudo scalar meson made 
up of one charm quark and another quark, for different (3 and continuum limit. For 
comparison, we also list the experimental values for the mesons in the last column . 



(3 


2.8 


2.6 


2.4 


continuum limit 


Experiment 


m(D) 


1835(35) 


1862(5) 


1860(9) 


1847(68) 


1869 


m(D s ) 


1976(21) 


1993(3) 


1959(5) 


2164(39) 


1969 


m(7] c ) 


2992(2) 


3037(1) 


3020(10) 


2843(23) 


2980 


m{D*) 


2015(27) 


1965(6) 


1993(33) 


2046(131) 


2007 


m(D*s) 


2140(18) 


2087(3) 


2080(13) 


2192(67) 


2112 



Finally, the results of the meson masses in the continuum are obtained by 
extrapolate their values at finite lattice spacings to the continuum limit. In 
Fig. 2, we show the mass ratios of various mesons over the mass of the J/ip 

2 The mass of J ftp (not included in Table 1) is by definition that of the experimental 
value since we have used the mass of J/i/j to determine the bare charm quark mass. 
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Fig. 1. Quadratic extrapolation/interpolation of My(m,m,u op t(m),i i pt(m)) versus 
the bare quark mass m at (3 = 2.6. The two curves, are the fits in the light and heavy 
quark mass regions, respectively. The corresponding fitting qualities and ranges are 
also shown in the figure. The fitting ranges are self-adjusted to yield minimum x 2 
per degree of freedom. The horizontal and vertical lines indicate the value for the 
physical <1> and J/ip meson and the fixed value of bare quark mass and m^' hy \ 

respectively. 

versus lattice spacing a s . Linear fits are performed for these mass ratios and 
the results are indicated by the straight lines. The extrapolated results are 
also indicated by the points at a s = 0. The numerical results are also listed 
in Table 1. The fitting quality for most cases are reasonable. However, since 
we have only 3 different values of /3 for this preliminary study, the continuum 
limit extrapolations suffer big errors. Needless to say, more points at different 
values of a s , preferably at smaller a s , will help to improve the final result at 
the continuum limit. 

To compare with the experiment, we have also tabulated the results from 
Particle Data Group (PDG) in the last column of Table 1. It is seen that 
the qualitative agreement is good, especially for D, D* and D* mesons. The 
mass values of the D s and r\ c are away from the experimental value by 5a. 
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Fig. 2. Continuum limit extrapolation for each meson mass over J /ip meson mass. 
Each panel indicates different meson. Data points (shown by green points) in the plot 
are ratios of the mass of each meson to the mass of J/ip for different (3 (corresponding 
to different value of a s /ro). The red points represent the physical mass ratios at the 
continuum limit. 

These disagreements might come from the fact that our lattice spacing is still 
somewhat large and the continuum extrapolation is not very well controlled. 
Of course, quenching effects are also present which might be important for the 
fine structure of the spectrum. 



3 conclusion 



In this paper, we present a preliminary study of the quenched meson spectra 
for the charmed mesons and charmonium states. To perform such a calculation 
using the tadpole improved anisotropic Wilson quark action, we tune the bare 
speed of light parameter v in a quark mass dependent way with quark mass 
values ranging from the strange to the charm. The optimal values of v are 
obtained for various values of (3 using the pseudo-scalar meson dispersion 
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relations. Using the tadpole improved anisotropic Wilson action with these 
optimized parameters, we determine the ground state meson spectra for D, 
D s , i] c , D* , D* and J/ip mesons. The results are extrapolated towards the 
continuum limit where qualitative agreement with the experiment is found. 
With the method outlined in this preliminary study, we plan to carry out 
spectrum calculations of other excited meson states, baryon states and possibly 
hybrid states involving a charm quark in the future. 
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